Electron spin diffusion at the interface of multiferroic oxides 
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We study the spin diffusion in a two-dimensional electron gas at the interface of oxide heterostruc- 
ture LaAIOs/SrTiOa grown on multiferroic TbMnOa at 15 K by means of the kinetic spin Bloch 
equation approach. The spiral magnetic moments of Mn^"*" in TbMnOs interact with the diffusing 
spins at the LaAlOa/SrTiOs interface via the Heisenberg exchange interaction. It is demonstrated 
that the spin diffusion length is always finite, despite the polarization direction of the injected spins. 
Our study also reveals the important role played by the Coulomb scattering, which can effectively 
suppress the spin diffusion. 
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I. INTRODUCTION 

During the past decades, spin-based electronics has 
been a growing area of research due to its promising 
apphcations.^^'^ Among various electron systems, the 
two-dimensional electron gas (2DEG) has been widely 
studied in spin injection, spin relaxation, spin transport 
and spin Hall effect, etc..^^^'^ Most of the 2DEG sys- 
tems are achieved in the semiconductor quantum wehs 
or heterostructures.^* The strict 2DEG, e.g., the one ex- 
isting in a graphene monolayer, also attracts much at- 
tention recently.^-'^"^'*^^^"^^ Besides, the 2DEG has also 
been experimentally realized at the interface of insu- 
lating oxides such as LaAlOa/SrTiOs,^^"^'* opening the 
way for the oxide-based nanoelectronics. Starting from 
this oxide heterostructure, very recently Jia and Be- 
rakdar proposed a trilayer system for new functionali- 
ties, which is grown as LaAlOa/SrTiOa/TbMnOa along 
the c-axis of orthorhombic TbMnOa. ^^"^^ As a multifer- 
roic material, TbMnOa can exhibit both magnetic and 
electric orders simultaneously and provide a unique op- 
portunity to exploit the multifunctionality of a single 
material.^^"'^^ SrTiOs in the middle of the trilayer sys- 
tem is a few layers thick so that the 2DEG at the in- 
terface of LaAlOs/SrTiOs can be affected by the mag- 
netic order in TbMnOs. Based on this structure, Jia and 
Berakdar predicted the appearance of a persistent spin 
current in the 2DEG due to the spiral geometry of the 
local magnetic order. They also proposed a flash mem- 
ory model based on the claim that the injected spins 
in the 2DEG have no decay of spin polarization along 
the a-axis when travelling along the 6-axis.^® This claim 
was inferred^^ from the fact that the spin polarization 
along the [110] direction does not decay in (001) GaAs 
quantum-well-based 2DEG with identical Rashba'^^ and 
Dresselhaus'^'^ spin-orbit coupling strengths, in the ab- 
sence of any external magnetic field. Nevertheless, 
this claim is questionable in the 2DEG formed at the in- 
terface of the oxide, because the local magnetic moments 
serve as a static magnetic field in the Voigt configuration 
in the spiral frame^^ and, as pointed out by Weng and 



Wu,^^ a magnetic field in the Voigt configuration can 
effectively suppress spin diffusion/transport even when 
the D'yakonov-Perel' spin-orbit coupling'^^ is absent. '^^''^^ 
In this work, we investigate the spin diffusion along the 
6-axis at the interface of LaAlOs/SrTiOs closely above 
TbMnOa by means of the kinetic spin Bloch equation 
(KSBE) approach.^ It is found that the spin diffusion 
length is finite in spite of the spin polarization direction 
and the Coulomb scattering makes marked effect on the 
spin relaxation in diffusion . 

This paper is organized as follows. In Sec. II, we give 
the Hamiltonian of the 2DEG in the collinear frame with 
parallel to the c-axis, and in the spiral frame with &z 
parallel to the local magnetic moment, respectively. We 
then present the KSBEs in the spiral frame. In Sec. Ill 
we study the spin diffusion based on the KSBEs. We 
summarize in Sec. IV. 



II. HAMILTONIAN AND KSBES 

The 2DEG formed at the LaAlOa/SrTiOa interface of 
LaAlOa/SrTiOa/TbMnOa compounds is schematically 
shown in Ref. 25. We assume it is confined by an in- 
finitely deep square potential well with width a = 3 nm.^^ 
The coordinate system is set as x = b, y = —a and z = c. 
When the temperature T is lower than the ferroelectric 
Curie temperature (^ 27 K), Mn3+ magnetic moments 
in TbMnOa form a spiral order with the local magnetic 
moment M(r) = sin(q • r)x -I- Mz cos(q • r)z, where 
q = (0.27 X X'O'T^) is the modulation vector.29.4i-43 
h = 0.586 nm and c = 0.749 nm are the lattice parameters 
of TbMnOa. 2^ In reality ^ (e.g., M^/M^ ~ 1.4 
at T = 15 K),^^'"'^"^'^ however for simplicity we take 
^ Mz = M following Jia and Berakdar. ^■'^'^^ It 
is assumed that the 2DEG only interacts with the lo- 
cal magnetic momentums of Mn'^"'" on the surface of 
TbMnOa (i-e., the plane z = 0) via the Heisenberg ex- 
change interaction. ^^'^^ As a result, the Hamiltonian of 
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the 2DEG veads^^^^^ 



p2 

H = h Jfir 

2m* 



(1) 



with the first term on the right-hand side of the equa- 
tion representing the kinetic energy and the second 
term, the exchange interaction. m* is the effective 
electron mass set to be lOme (mg is the free-electron 
mass),^^'^^ J stands for the coupling strength, fir = 
{sm{qxx),0,cos{qxx)) denotes the unit vector along the 
local magnetic moment located on ^; = plane, P is 
the momentum operater and cr is the vector of the Pauli 
matrices. Performing a local rotation around the y-axis 
in the spin space as .H" = U^{x)HUg{x) with Ug{x) = 
^-iq^xdyji^ one obtains in the spiral frame^^'^^ 

This Hamiltonian can be written in the momentum space 

as 



in which 



(3) 



(4) 



^2 2 

and a uniform energy displacement eo = is omitted. 
The position-dependent 3x3 orthogonal rotation matrix 
R{x)-, which obeys U^{x)cTUg{x) = R{x)ct, reads 



Rix) = 



cos{qxx) sin(g'a,x) 

10 
— sm{qxx) cos(g'xa;) 



(5) 



R'^{x) transforms any spin- vector in the spiral frame 

back to the collinear frame. Our study is performed by 
first solving the KSBEs in the spiral frame and then ob- 
taining the spin diffusion properties in the collinear frame 
with the aid of R?'{x). 

In the study, polarized spins are injected at x = and 
diffuse along the .x-axis. The system is uniform along 
the fixed j/-axis in both frames. In the spiral frame, the 
KSBEs read^ 



dpi^{x, t) _^ e 9^(a;,t) dpkix,t) 
dt H dx dkx 

-I- - [hk • <T, pu{x, t)\ + — 



m" 



dx 



scat 



= 0, 



(6) 



where pi^{x,t) are the single-particle density matrices of 
electrons with wave-vector k at position x and time t. 

f ^'^irc'*'' ^^'^ ^^'^ driving terms with the electric 

potential satisfying the Poisson equation V^\E'(x,t) = 
e[Ne{x,t)-No\/ (aKoeo), where Ne{x,t) = J2k^[pk{x,t)] 
stands for the electron density at position x and time t, 
No is the background positive charge density and kq 24 



is the relative static dielectric constant. In the ab- 
sence of the external electric field, the driving term is 

negligible. p^{x,i)] and ^ ^^^^ 



in Eq. (6) are the diffusion, coherent and scattering 
terms, respectively. We lack detailed information of the 
electron-phonon and electron-magnon scatterings in this 
material, except that the electron-phonon scattering in 
SrTiOa was suggested to be typically very weak.^'* How- 
ever, as we are interested only in the low temperature be- 
havior (we take T = 15 K< Tc in this paper), both scat- 
terings are neglected. The electron-impurity scattering 
is also neglected in the high-mobility 2DEG (for electron 
density A'^e 10"^^ cm~^ and mobility ^ 10^ cm^/(V-s).^^ 
the impurity density is calculated to be ~ 10* cm~-^, cor- 
responding to which the electron-impurity scattering is 
found to be negligible). The scattering left in our con- 
sideration is the electron-electron Coulomb scattering, 
which has been proved to be important in spin relaxation 
and spin transport in semiconductors. ^'^^""^^ The detailed 
expression of the Coulomb scattering can be found in 
Ref. 49. By solving the KSBEs one obtains the spin dif- 
fusion properties from pk(a;, +oo), the steady-state dis- 
tribution of density matrices. 



III. SPIN DIFFUSION 
A. Spin relaxation in time domain 

Equation (4) is analogous to the spin-orbit coupling 

in (001) GaAs quantum- well-based 2DEG with identi- 
cal Rashba and Dresselhaus spin-orbit coupling strengths 
(only the linear Dresselhaus term is considered) when J is 
set to zero. There, when the x-axis is set along [110], the 
spin-orbit coupling reads hk = {0,2"fkx,0), with 7 be- 
ing the coefficient of the Rashba (or Dresselhaus) spin- 
orbit coupling. ^^'^^'^*'^° Therefore, if the spin polariza- 
tion is along the y-axis, there is no spin precession and 
hence the spin relaxation is suppressed. ■'■''•■"^■"'*'''" How- 
ever, the exchange interaction in oxide provides a static 
magnetic field J along the 2-axis in the spiral frame, as 
shown in Eq. (4) . This field rotates the spin polarization 
away from the j/-axis. The spin polarization away from 
the 2/-axis feels the inhomogeneous broadening, i.e., 
the momentum-dependent spin-orbit coupling from the 
y-component of hk, which, together with any scattering, 
leads to irreversible spin relaxation. ^'''^'^^ Moreover, the 
effective magnetic field mixes spin relaxations along the 
x-z plane to the one along the y-axis, similar to the case 
based on semiconductor quantum wells. Therefore, 
the spin relaxation time along the y-axis is not infinite 
as claimed by Jia and Berakdar,^^ neither does the spin 
relaxation along other directions. The above scenario 
also happens to (110) GaAs quantum wells. The effec- 
tive magnetic field induced by the Dresselhaus spin-orbit 
coupling is oriented along the quantum-well growth di- 
rection when only the lowest subband is relevant, which 
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leads to an infinite spin relaxation time when the spin 
polarization is along that direction. However, by apply- 
ing a static magnetic field in the quantum- well plane, the 
spin relaxation time becomes finite. ^^'^^ 

To be quantitative, we numerically solve the KSBEs 
in time domain^'^*'"'^"'*^ with the Coulomb scattering ex- 
plicitly included. In Fig. 1, the time evolution of spin 
polarization (with an initial value 5 %) is plotted. It is 
clearly shown that when the initial spin polarization n is 
along y, the spin polarization relaxes when J ^ 0. How- 
ever, for other spin polarization directions such as z, spin 
relaxes even when J = 0. 
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FIG. 1: (Color online) Time evolution of spin polarization 
along different directions under different J (in the figure J = 



B. Spin diffusion: analytical and numerical study 

In this section we investigate the spin diffusion of the 
2DEG. We first consider a simplified case without the 
scattering for which the KSBEs can be solved analyti- 
cally. Then we numerically solve the KSBEs in the pres- 
ence of the Coulomb scattering. 

When the scattering terms are excluded (and the driv- 
ing terms are absent), the KSBEs can be solved in the 
steady state as 



_ -iUJk-<Tx/2 



with 



2m* 



-hk 



/9k(0, +(X))e 



2m* , 



iUJk-Crx/2 



(7) 



(8) 



depicting the spin precession frequency in spatial domain. 
It has been shown that the inhomogeneous broadening in 
spin diffusion/transport is determined by Wk rather than 
hk in the time domain, ^■'^^'^^''''''^^ with arising from 
the diffusion term. It becomes evident that J, playing 
the role of a static magnetic field in the time domain, 
now leads to the inhomogeneous broadening due to 
and hence a finite spin injection length. Similar effect 
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FIG. 2: (Color online) S'S (a;)/S'S (0) vs. x for different fi and 
J. (a) and (d): n = x; (b) and (e): n = y; (c) and (f): n = z. 
(a)-(c) are plotted in a large scale of x while (d)-(f) are in a 
small scale. 



was first predicted theoretically by Weng and Wu back 
in 2002 (Ref. 38) and realized experimentally in bulk 
silicon. Moreover, due to the presence oiqx in Eq. (8), 
the spin diffusion length is finite even when the injected 
spins are polarized along the z-axis. 

The spatial distribution of spin polarization is obtained 
from the density matrices given by Eq. (7). We consider 
the states with fc^; > because only these states can prop- 
agate along the a;-axis and the states with fc^: < are not 
spin-polarized in the absence of scattering.^® Therefore 
the density matrices with > at x = are fixed 
as boundary conditions. The spin polarization Sk(a:^) = 
Tr[pk(2;, -)-oo)<t] is calculated to be Sk(a:^) — Fk(a;)Sk(0), 
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with 



ft2 



/ik,z sin(wka;) 
cos(ti;k. 
cos(ajka;)] 



hia^y sin(wka;) 
^k,y^k,z[l - cos(a;ka;)] 



(9) 



in which uj^ = |u;k|, /Jk,!y(z) = /''k.y(2)/|hk| and Sk(0) = 
Tr[pk(0, 0)cr]. The total spin signal in the spiral frame is 

S(x) = Sk(-^) = E ^k(a:)Sk(0), (10) 

fcx>0 fcx>0 

and in the collinear frame reads 

S"(a;) = R'^{x)S{x). (11) 
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FIG. 3: (Color online) S^{x) / S^(0) vs. x for different n when 
Qx is set to be zero. J = 0.2. 

In Fig. 2, the projection of spin signal S''{x) on the 
injected spin polarization direction fi, S^{x) = S'^lx) ■ n, 
calculated from Eq. (11) is plotted against x with differ- 
ent J [the data are scaled by S'S(O)]. In the computa- 
tion, we set Sk(0) = [/(ek - M+) - /(^k - /i-)]n.37.4o,48 
/(^k — M{) with Ek — h^k^/2m* is the Fermi distribution 
under temperature T. The chemical potential /i^ is de- 
termined by Ne, the electron density of the 2DEG, and 
Ps, the polarization of injected spins, at a; = via the 
relation J2k /(^k - A^c) = 0.5(1 -I- ^Ps)Ne.^'^^'^°'^^ We set 
Ne = 10" cm-2 (Ref. 22) and ^^=5 %. fi is set as x, 
y and z. It is shown that when J 7^ 0, the spin polar- 
ization decays during spin diffusion. When J increases, 
the decay becomes faster. It is noted that Si{x) relaxes 
to a finite value for n = y in the absence of scattering. 
However, when the scattering is included, the spin polar- 
ization will completely relax to zero during spin diffusion 
(this will be shown later). Moreover, the spatial distribu- 
tion of spin polarization shows high-frequency oscillations 



when J 0, as detailed in Figs. 2(d)-(f) in a smaller scale 
of X. The frequency of these oscillations is comparable 
to the rotation frequency of the spiral magnetic moments 
and insensitive to the exchange coupling strength. For 
comparison, we recalculate Si^{x) by replacing the spiral 
order with the ferromagnetic order, i.e., aligning all the 
magnetic moments uniformly along the z-axis. Hence 
qx = 0. The corresponding results are shown in Fig. 3, 
where the high-frequency oscillations disappear and the 
spin polarization along the z-axis does not decay due to 
the absence of the inhomogeneous broadening. 
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FIG. 4: (Color online) S^{x)/S^{0) vs. x for different fi and 
J. Here Sf^{x) are obtained by numerically solving the KSBEs 
in the presence of the Coulomb scattering. 

It is noted that the above decay of spin polarization 
along the spin diffusion in the absence of scattering is 
due to the interference effect.'^* Now we numerically solve 
the full KSBEs with the Coulomb scattering explicitly in- 
cluded, facilitated with the double-side injection bound- 
ary conditions. For the cases with J 7^ 0, we plot 
Si as a function of position x in Fig. 4. It is shown that 
the Coulomb scattering effectively leads to the spin relax- 
ation along the spin diffusion. By comparing Figs. 4 and 
2, it is observed that with the scattering, the spin diffu- 
sion is suppressed along all directions.'*" Especially, the 
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spin polarization along the y-axis relaxes to zero com- 
pletely during the spin diffusion. It is also noted that 
the amplitude of the high-frequency oscillations shown 
in Fig. 2 becomes much smaller. 

IV. CONCLUSION 

In conclusion, we have studied the spin diffusion 
in a 2DEG at the interface of oxide heterostructure 
LaAlOa/SrTiOs closely above the multiferroic TbMnOs 
by means of the KSBE approach. At low temperature 
(< 27 K), Mn"^+ magnetic moments in TbMnOs form a 
spiral order and couple with the 2DEG via the Heisenberg 
exchange interaction. It is found that in the presence of 



the spiral magnetic order, in contrast to the claims in the 
literature, the spin polarization decays during spin diffu- 
sion, despite the polarization direction of injected spins. 
We also report that the electron-electron Coulomb scat- 
tering suppresses spin diffusion effectively at low temper- 
ature. 
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